Suppression of Transport of an Interacting Elongated Bose-Einstein Condensate in a 

Random Potential 
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We observe the suppression of the ID transport of an interacting elongated Bose-Einstein conden- 
sate in a random potential with a standard deviation small compared to the typical energy per atom, 
dominated by the interaction energy. Numerical solutions of the Gross-Pitaevskii equation repro- 
duce well our observations. We propose a scenario for disorder- induced trapping of the condensate 
in agreement with our observations. 
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Atomic Bose-Einstein condensates (BEC) in optical 
potentials are a remarkable system in which to revisit 
standard problems of condensed matter physics, e.g. su- 
perfluidity and quantum vortices, the superfluid to Mott 
insulator transition, or Josephson arrays 0. Another 
important topic in condensed matter physics is that of 
transport in disordered materials, with relevance to nor- 
mal metallic conduction, superconductivity and super- 
fluid flow in low temperature quantum liquids. This is 
a difficult problem and it has led to the introduction of 
intriguing and non-intuitive concepts, e.g. Anderson lo- 
calization 0, |s| , percolation Q and Bose and spin 
glasses. It also has a counterpart in wave physics, e.g. 
in optics and acoustics, specifically coherent diffusion in 
random media 0. The main difficulty in understanding 
quantum transport arises from the subtle interplay of in- 
terference, scattering onto the potential landscape, and 
(whenever present) interparticle interactions. 

Transport properties of BECs in periodic optical lat- 
tices have been widely investigated, showing lattice- 
induced reduction of mobility HlH, 0] and self-trapping 
Within the context of random potentials, most of 
the recent theoretical efforts 12] have considered disor- 
dered or quasi-disordered optical lattices where a large 
variety of phenomena have been discussed such as the 
Bose-glass phase transition |F3j , localization [F3J, \VA\ , and 
the formation of Fermi-glass, quantum percolating and 
spin-glass phases in Fermi-Bose mixtures 

[H E3 Ef- 
fects of disorder on BECs have also been addressed in 
connection to superfluid flows in liquid helium in porous 
media In particular, the depletions of the conden- 
sate and of the superfluid fractions have been calculated 
in Ref. |l7| , and a significant shift and damping of sound 
waves have been predicted in Ref. [l^. Apart from the 
(undesired) fragmentation effect of a rough potential on 
trapped cold atoms and BECs on atom chips [lflj . there 
are few experiments on BECs in random potentials |2(ij . 



In this Letter we report on the strong reduction of mo- 
bility of atoms in an elongated BEC in a random poten- 
tial |2l|. Starting from a BEC in a 3D highly elongated 
harmonic trap, we turn off the axial trapping potential 
while maintaining strong transverse confinement, and we 
monitor both (i) the axial expansion driven by the repul- 
sive interactions and (ii) the motion of the center of mass 
of the BEC. When the BEC is subjected to a ID random 
potential created by laser speckle, the axial expansion 
is strongly inhibited and the BEC eventually stops ex- 
panding (see Fig. The final rms size L decreases as 
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Figure 1: Time evolution of the axial rms size L of the BEC, 
for various amplitudes cry of the random potential, all smaller 
than the chemical potential \x [7 = cry// 1 — (o), 0.2 (•) and 
0.7 (■)]. The axial trapping frequency is initially u z /2n = 
6.7 Hz and is relaxed during the first 30 ms (lu z t < 1.26) of 
the expansion time (grey band). Each point corresponds to 
an average over three measurements; error bars represent one 
standard deviation. The solid lines are linear fits to the data 
and the dashed lines are guides to the eye. Inset: Motion of 
the center of mass of the BEC during axial expansion for the 
same values of 7. Both sets of data show a strong suppression 
of transport of the BEC in the presence of disorder. 
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Figure 2: Optical setup used to create the random speckle 
potential. The BEC is at the focus of the lens system with 
its long axis oriented along the z direction. 

the standard deviation ay of the random potential in- 
creases. The same effect has been observed for various 
realizations of the random potential. We also observe 
that the center of mass motion provoked by a longitudinal 
magnetic 'kick' at the time of release is strongly damped 
and is stopped in about the same time (see Fig.QJ . These 
observations are not made in a regime of tight binding, 
i.e. we observe this localization effect [3lJ for amplitudes 
of the random potential small compared to the chemi- 
cal potential. One may wonder whether our observations 
can be interpreted in terms of Anderson localization Q • 
In fact, in our situation, the interaction energy plays a 
crucial role, and the healing length is smaller than the 
typical distance between the speckle grains. This implies 
a different scenario, which we discuss in this Letter. 

We create an elongated 87 Rb BEC in an iron-core elec- 
tromagnet Ioffe-Pritchard trap jl^, with oscillation 
frequencies, oj±/2it = 660(4) Hz radially and uj z /2tt = 
6.70(7) Hz axially. BECs of typically 3.5 x 10 5 atoms 
are obtained, with Thomas-Fermi half-length Ltf = 
150 /im and radius Rt v = 1.5 /mi, and chemical poten- 
tial /j,/2ttH ~ 5 kHz 32]. The random potential is turned 
on at the end of the evaporative cooling ramp and we 
further evaporate during 200 ms to ensure that the BEC 
is in equilibrium in the combined harmonic plus random 
potential at the end of the sequence. 

To create the random potential, a P < 150 mW blue 
detuned laser beam with optical wavelength A ~ 780 nm, 
perpendicular to axis z is shone through a scattering 
plate and projects a speckle pattern [2|| on the BEC (see 
Fig. [21 • The scattered beam diverges to an rms radius of 
1.83 mm at the BEC. 

A speckle field is defined by (i) a random intensity 
J(r) with exponential statistical distribution for which 
the standard deviation equals the average intensity 07 = 
(J) and (ii) an intensity correlation length Az, defined as 
the 'half- width' of the autocorrelation function [2(| : 

Az = 1.22 Xl/D, (1) 

where D is the beam diameter at the scattering plate 
and I is the distance from the lens to the BEC. We ob- 
serve the speckle intensity distribution on a CCD cam- 
era placed at the same distance as the atoms. From 



this, we determine the autocorrelation function to obtain 
the grain size Az for various beam diameters D. Tak- 
ing into account the modulation transfer function of 
the camera, we find that the measured grain size obeys 
Eq. to within 2%. For our setup (I = 140(5) mm and 
D = 25.4(1) mm), Eq. © gives Az = 5.2(2) fjm. This 
is an order of magnitude greater than the healing length 
£ = (8nna)~ 1/2 = 0.11 fim of the trapped BEC. Since 
i?TF < Az <C Ltf, the optical potential is effectively ID, 
with the trapped BEC spread over about 45 — 50 wells 
in the axial direction. We characterize the amplitude of 
the random potential ay with respect to the chemical 
potential fi by |2sj : 

_ ay _ 2 T 2 ai 1 _ 1 T 2 ai / 15aN \ ~ 2/5 
H 3 26 I s /x w6i/ s \ ai^ / 

with U=(u> z ui'j_) 1 / 3 and a] 10 =(h/muj) 1 ^ 2 , m the atomic 
mass, N the BEC atom number, Is — 16.56 W/m 2 the 
saturation intensity, r/27r = 6.01MHz the linewidth, 
a = 5.31 nm the scattering length and 6 the laser de- 
tuning (between 0.15 nm and 0.39 nm in wavelength). 
The factor 2/3 accounts for the transition strength for 
7r-polarized light. Taking into account our calibration 
uncertainty, we measure 7 within ±20%. For our param- 
eters, the spontaneous scattering time l/r sc is always 
larger than Is, i.e. much longer than the experiment. 

To study the coherent transport of the BEC in the 
random potential, we open the axial magnetic trap while 
keeping the transverse confinement and the random po- 
tential unchanged. After lowering the current in the ax- 
ial excitation coils, the axial trapping frequency oj z /2tt 
is smaller than 1 Hz |33|, Opening the trap abruptly in- 
duces atom loss and heating, therefore the trap is opened 
in 30 ms to avoid these processes. Once the current in the 
axial coils has reached its final value we have a BEC of 
N ~ 2.5 x 10 5 — 3 x 10 5 atoms in the magnetic guide. 

After a total axial expansion time r (which includes the 
30 ms opening time) , we turn off all remaining fields (in- 
cluding the random potential) and wait a further 15 ms of 
free fall before imaging the atoms by absorption. During 
this time-of- flight, the axial rms size of the BEC does not 
increase more than 5%. From profiles of the absorption 
images we evaluate the axial rms size L [34| which we plot 
in Fig.^versus the axial expansion time r. In the absence 
of the random potential (7 = 0), we observe that the rms 
size L grows linearly at a rate wrms ~ 2.47(3) mms -1 in 
agreement with the scaling theory |2^|. In the presence 
of the random potential, the expansion dynamics changes 
dramatically. For a sufficiently high amplitude, the ex- 
pansion is significantly reduced and the BEC eventually 
stops expanding. In addition, we observe the damping 
of longitudinal motion of the center of mass of the BEC 
(see inset of Fig.GJ. This motion is triggered by an axial 
magnetic 'kick' during the opening of the trap. 

These results show a transition from non-inhibited to 
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Figure 3: Rms size L of the BEC versus 7 after an axial 
expansion time uj z t = 4.84 (r = 115ms). The open circles 
correspond to the curves of Fig. Q 

inhibited transport as the speckle amplitude is increased. 
This is studied in further detail by measuring the BEC 
rms size after a fixed axial expansion time of 115 ms 
(uj z t — 4.84) for different amplitudes ay of the random 
potential. The results are shown in Fig. [3] We see that 
above a value 7 = 0.15, the rms size decreases with 7. 

From the absorption images, we also evaluate the den- 
sity in the magnetic guide, after correcting for radial ex- 
pansion during the time-of- flight. We observe that the 
density at the center of the BEC does not drop by more 
than a factor of 2 for 7 > 0.2. Therefore, we conclude 
that the interaction energy dominates at the center of the 
BEC trapped by disorder, a point we shall discuss below. 

To understand the suppression of expansion of the 
BEC in the random potential, we have performed nu- 
merical calculations of the BEC dynamics in the mean- 
field Gross-Pitaevskii approach. We consider a BEC 
trapped in a cylindrically-symmetric 3D-harmonic trap 
with frequencies lu±_ and uj z in the radial and axial di- 
rections respectively. Assuming tight radial confine- 
ment (hu>± 3> huj z , fi, fceT), the dynamics is reduced 
to ID. In addition, the BEC is subjected to a static 
random potential V(z) = ovv(z) where v(z) is a nor- 
malized numerically-generated speckle pattern j2|| with 
(v) 2 = (v 2 )/2 = 1. This slightly differs from the exper- 
imental situation where the BEC is very elongated but 
not strictly ID. However, in the experiment, the BEC 
is guided in a ID random potential so that the radial 
size only slightly changes and, due to the different time 
scales in the axial (1/oj z ) and radial (1/w_l <C l/w 2 ) di- 
rections, the radial size adapts adiabatically to the ax- 
ial size. Thus we expect that the ID simplified model 
captures the physics of the experiment. We consider pa- 
rameters close to the experimental situation (see above). 
In particular, the healing length (£ ~ 8 x 10~ 4 Ltf) an d 
the speckle correlation length (Az ~ 0.049Ltf) are much 
smaller than the size of the BEC. 

We first compute the static ID BEC wavefunction in 
the combined (harmonic plus random) trap. Because 
£ <C Az, the density profile simply follows the modula- 
tions of the combined trap in the Thomas-Fermi regime: 
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Figure 4: a) Time evolution of the rms size L of the BEC in 
the random potential V(z) for various values of the speckle 
amplitude av = 7/i as obtained from the numerical cal- 
culations, b) Density profile (black) and random potential 
V{z)/g\T> (gray) for 7 = 0.2 at ui z t = 10. c-d) Enlargement 
of density profile at uj z t — 10 (black solid) and uj z t = 20 
(blue dotted). The red dashed line in c) is the Thomas- Fermi 
prediction (see text). 

I 2 = [A* — moj 2 z 2 /2 — V(z)]/giD in the region where 
fi > muj 2 z 2 /2 + V(z) and | -0 (2;) | 2 = elsewhere. Here, 
to is the atomic mass and giu — 2hauj± the ID interac- 
tion parameter. At time r = 0, we suddenly switch off 
the axial harmonic confinement while keeping unchanged 
the interaction parameter giD and the random potential 
and we compute the time-evolution of the BEC. The re- 
sults for the axial rms size L of the BEC are plotted in 
Fig. ^Jl for various amplitudes of the random potential. 
In the absence of disorder, the evolution of the BEC cor- 
responds to self-similar expansion with scaling parameter 
b(t) ~ V2uj z t 0. In the presence of disorder (7 > 0.15), 
after initial expansion, the BEC stops expanding. This 
is qualitatively the same behavior we observed in the ex- 
periment. The quantitative agreement is also reasonably 
good. For example, for 7 = 0.2, the BEC expands by a 
factor of ~ 4 in the numerics (~ 3 in the experiment) and 
is trapped after a transient expansion time of lo z t ~ 8 
(uj z t ~ 6). This strong suppression of expansion corre- 
sponds to disorder-induced trapping of the BEC. 

We now describe a scenario for disorder-induced trap- 
ping of the BEC after release of the axial harmonic con- 
finement. For small enough amplitudes of the random 
potential, the initial stage of expansion can be described 
using the scaling theory j^. According to this, the fast 
atoms populate the wings of the expanding BEC whereas 
the slow atoms are close to the center. It is thus tempt- 
ing to distinguish two regions of the BEC: (i) the cen- 
ter where the interaction energy dominates the kinetic 
energy and trapping is due to the competition between 
interactions and disorder, and (ii) the wings where the 
kinetic energy exceeds the interactions and trapping is 
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rather due to the competition between the kinetic energy 
and disorder. 

In the center, the average density and thus the effective 
chemical potential ~p slowly decrease during the expan- 
sion stage. As the interaction energy is much larger than 
the kinetic energy, the local density adiabatically follows 
the instantaneous value of ~p in the Thomas-Fermi regime: 
|?/>(z)| 2 = [71 — V(z)]/gm in the region where ~p > V(z) 
and |^(z)| 2 = elsewhere. This agrees with our nu- 
merical results (see Fig. St). This evolution stops with 
fragmentation, i.e. when the BEC meets two peaks of the 
random potential with amplitudes larger than ~p. Using 
the statistical properties of the random potential [2(J , we 
can estimate the probability of such large peaks and we 
conclude that this happens when the central density no 
reaches the value 



n ~ 1.25 ( — ) In 
gin 



0.47L T f 
Az 



(3) 



This formula is in good agreement with our numerical 
and experimental findings 30]. 

Due to the small density, the situation is completely 
different in the wings which are populated by almost 
free particles interacting with the disordered potential. 
The BEC wavefunction thus undergoes disorder-induced 
multiple reflections and transmissions and is ultimately 
blocked by a large peak of the speckle potential. There- 
fore, the BEC is not in the Thomas-Fermi regime and 
the local density is not stationary (see Fig. 0Ji). Due to 
conservation of energy, the kinetic energy per particle e 
is of the order of the typical energy in the initial BEC 
(e ~ jti) so that the typical wavelength of the fluctuations 
in the wings is of the order of the healing length in the 
initial BEC A w ~ £ = h/i/2mji. 

This scenario of disorder-induced trapping is accu- 
rately supported by our numerical integration of the 
Gross-Pitaevskii equation. In particular, the density pro- 
files plotted in Fig. 0] show the static Thomas-Fermi 
shape in the center and time-dependent fluctuations in 
the wings with typical wavelength A w ~ £ [3(J ■ 

In conclusion, we have experimentally investigated 
transport properties of an interacting BEC in a random 
potential. Controlling the strength of disorder, we have 
observed the transition from free expansion to absence 
of diffusion as disorder increases. We have presented 
numerical simulations that reproduce well the observed 
suppression of expansion and we have discussed a the- 
oretical model that describes the scenario for disorder- 
induced trapping. In the future, it would be interesting 
to further investigate this highly controllable system, for 
example by changing the correlation length of disorder or 
employing Bragg spectroscopy to probe the momentum 
spectrum of the BEC 0. 
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